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ABSTRACT 

We assess the detection prospects of a gravitational wave background associated with 
sub-luminous gamma-ray bursts (SL-GRBs). We assume that the central engines of a 
significant proportion of these bursts are provided by newly born magnetars and con- 
sider two plausible GW emission mechanisms. Firstly, the deformation-induced triaxial 
GW emission from a newly born magnetar. Secondly, the onset of a secular bar-mode 
instability, associated with the long lived plateau ob served in the X-ray afterglows 
of many gamma-ray bursts (|Corsi fc Me szaros 2009a). With regards to detectability, 
we find that the onset of a secular instability is the most optimistic scenario: under 
the hypothesis that SL-GRBs associated with secularly unstable magnetars occur at 
a rate of (48 — 80) Gpc _3 yr _1 or greater, cross-correlation of data from two Einstein 
Telescopes (ETs) could detect the GW background associated to this signal with a 
signal-to-noise ratio of 3 or greater after 1 year of observation. Assuming neutron 
star spindown results purely from triaxial GW emissions, we find that rates of around 
(130 — 350) Gpc yr _1 will be required by ET to detect the resulting GW background. 
We show that a background signal from secular instabilities could potentially mask 
a primordial GW background signal in the frequency range where ET is most sen- 
sitive. Finally, we show how accounting for cosmic metallicity evolution can increase 
the predicted signal-to-noise ratio for background signals associated with SL-GRBs. 

Key words: gravitational waves - gamma-ray bursts - supernovae: general - cos- 
mology: miscellaneous 



1 INTRODUCTION 

The two closest recorded gamma-ray bursts (GRBs), GRB 
980425 (36 Mpc) and GRB 060218 (145 Mpc), along with 
GRB 031203 , asso ciated with a host galaxy at ~ 480 Mpc 
l|Feng fc Fcodl2010t ). make up a sub-class of long duration 
GRb£\ (LGRBs) known as sub-luminous GRBs (SL-GRBs) 
JCobbet al.ll2006l: iLiang et al.ll2007l : iGuetta fc Delia Valid 
120071 : IVirgili et al.ll2009h . This class of GRB have isotropic 
equivalent 7-ray energy emissions typically several or- 
ders of magnitude be l ow th o se of standard long-duration 
GRBs ilMurase et all 120061 : IGuetta fc Delia Valid 120071 : 



llmerito et al. 20081 ) suggesting that they could form 



unique population of bursts. 



* E-mail:email@address 

1 Hereafter, we refer to long GRBs as those having a T90 du- 
ration > 2s in agreement with the traditional classification by 



Observations have confirmed that at least some 
LGRBs are associated w ith the deaths of massive stars 
|Wooslev fc Bloomi r2006h One scenario for the g eneration 
of LGRBs is described by the collapsar model (|Wooslevl 
Il993l : iMacFadven et afll2o'oil h In this model, the inner part 
of a progenitor star (a Wolf-Rayet, WR, star) collapses, via 
a Type Ib/c supernova, forming a rapidly rotating black 
hole. High angular momentum enables the infalling matter 
to form an accretion disk, which in turn powers an ultra- 
relativisti c jet (a "fireball") that blasts through the stellar 
envel ope (jMeszaros fc Reesill992l : IWooslevlll993l ; ISari et all 
Il998h . A number of authors have suggested that, at least 
in some cases, GRB explosions may end in the formation 
of a highly magnetized n eutron star (NS), i.e. a magnetar 
rather than a black hole dUsovlll992l; Duncan fc Thompson! 
19921: iDai fc Lull 19981: lNabdl2007flBucciantini et all 12009} 



Zhang & Dai 2009). Additional support for this scenario has 



ratio n > 2s m agrcemeni 
iKouveliotou et all lll993l l 



come from detailed modeling of the spectra and light curve 
of SN 2006aj (Type Ic), associated with GRB 060218. This 
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analysis suggested that the explosion energy and ejected 
mass originated from a progenitor star with a zero-age main 
sequence mass of ~ 2 0Ma , implying the bi r th of a NS rather 
than a black hole (|Mazzali et alj 120061 ; ISoderberg et al.l 
120061 ). The fact that this burst was of the SL-GRB class 
suggests that a pro portion of such bursts may well be po w- 
ered by magnetars (|Toma et al.ll2007l ; iMurase et al.ll2006n . 

The magnetar scenario for GRBs has been invoked to 
explain recent observations by the Swift satellite Q, showing 
that a significant fraction of the X-ray afterglo ws of GRBs 
exhibit a shallow decay phase lasting 10 2 -10 4 s |zhanget al.1 
20061; iFan fc Xu 2006; Nousck et al. 20061; iLiang et al.ll2007l ; 



Yu et al.1 120071 ; lYama zaki 2009(). A number of stud ies have 



suggested that the long duration afterglow plateau may 
be powered by a newly born millisecond magnetar. This 
could channel slowly decreasing rotational energ y into a rel- 
ativistic outflow via magnetic dipole emis s ion dUsovlll992l; 



tag] 

Zhang fc Meszaros 2001; iFan fc Xul 120061; [PaH'Osso etall 



20ld ; lLvons et al.ll2010l ; [Yu et al.ll201(ih 



From the perspective of GW detection, the relative lo- 
cal proximity of observed SL-GRBs makes these sources 
appealing, and raises the possibility of mult i- messenger 
obse rvations by second and third generation GW detec- 
tors jKochanek fc Piranlll993l;lKobavashi fc Meszarosll2003l; 



Abbott et al.l |200S|; lAndersson et al.1 120091; iBloom et al.1 



20091 ; ICorsi fc Meszaroj l2009al lbl; lAbbott et al.l l2010bl lal) 

Any detection scenario, of course, depends on the frequency 
of events in the ne arby Universe. R a dio ob servations of 68 
local SNe Ib/c by ISoderberg et~aH {2006) show that less 
than 10% were associated with LG RBs. Based on the lo- 
cal rate of Type lb/ c sup ernovae l|Soderberg et alj 120061 : 
iGuetta fc Delia Valid [20071 ) this yields an extreme upper 
limit of ~ 2 x 10 3 Gpc" 3 yr _1 for SL-GRBs. Studies of SL- 
GRB rates over the last four years (to be discussed in the 
next section) have yielded estimates extending over a range 
(40 - 1800) Gpc" 3 yr _1 . These estimates are orders of mag- 
nitude greater than those of classical LGRBs. 

Although these rate estimates are encouraging, GW de- 
tections will also depend on the strength of the emissions. As 
LGRBs require rapid r otation to produce an accretion disc 
ilYoon fc Langerjl2005l ; IWooslev fc Jankal 120051 ; Ijanka et al.1 
120071 ; Hartman nl l2010l ). it is logical to consider post-collapse 
GW emission mechanisms equally dependent on rotation. 

A first possibility is that the strong magnetic fields 
of newly born magnetars , of order 10 14 — 10 16 G 
jDuncan fc Thompson! Il992l) . could lead to deformations 
that would dominate any flattening due to a fast rotatio n 
jOstriker fc Gunnll 19691 ; |Palomba|[200ol ; iKonno et~aill2000l ). 
If the deformation axis is offset from the spin axis, this 



could lead to GW emissions (Palomba 200 ll; ICutlerl 


2002; 


Stella et al. 2005 


; Regimbau & de Freitas Pachecol 


2006: 


DallOsso & Stella 


120071: iDallOsso et al. 20091). 




A second 


emission possibility, suggested 


by 



Corsi fc Meszaros! (|2009al lbD. is that GW emissions 
could accompany the electromagnetic dipole emissions 
of a newl y formed magnetar via a secular bar mode 
instability (IChandrasekharjll970| ; [Friedman fc Schutzlll978l; 



Lai fc Shapirol Il995l; lOu et al.l |2004| ; IShibata fc Karinol 



20041 ). As this instability occurs on a long timescale ~ 10 



10 4 s, it corresponds well with the observed X-ray plateau 
of some LGRBs. We note t hat we do not consider the GW 
emissi ons from R modes (|Owen et al. I QUI; iFerrari et al.1 
Il999bl ). as the effect of the magnetic field could suppress 
this GW ins tability in magnetars (|Rezzolla et al.1 l200ll ; 
lMendell|[200ll ). 

Despite the observed sample of SL-GRBs being small, 
and the origin of this class of bursts still some way from 
being clearly understood, the rate estimates based on cur- 
rent observations are of an equivalent order to those of 
other sources of pot entially detectable GW backgrounds , 
e.g. NS/NS mergers ([Regimbau fc de Freitas Pachecol 120061 ; 
lRegimbau1l2007l ; lRegimbau fc Mandidl200Sl ). For this reason 
we are motivated in this paper to determine if a GW back- 
ground from SL-GRBs, based on the two mechanisms out- 
lined above, could produce a detectable signal for advanced 
GW interferometric detectors such as ALIGC0 (Advanced 
Laser Interferometer Gravitational-wave Observatory) and 
VIRGC0 or third generation instruments such as the Ein- 
stein Telescope0 (ET). Although we expect there to be varia- 
tion in the single-source emission mechanisms, we take these 
models to represent average values. This assumption is rea- 
sonable, based on the fact that detection of a stochastic GW 
background can only yiel d information on the mean event 
emission of a population dRegimbau fc Mandid [2008). Ad- 
ditionally, we allow for uncertainties in both the event rates 
of SL-GRBs and in the frequency of occurrence of the two 
GW emission mechanisms considered by employing widely 
separated upper and lower limits. 

The possibility of a detectable continuous astrophysi- 
cal background signal is important, as it could mask the 
relic GW background signal from the earliest epochs of 
the Universe. Primordial backgrounds are expected to be 
produced by l arge numbers of dynamical events in the 
early Universe ijGrishchukl 1 19751 ) - this signal is expected 
to be isotropic, stationary and unpolarized. An upper limit 
on the energy density of the primordial GW background 
normalized by the critical energy density of the universe, 
was recently set as < 6.9 x 10 -6 in the frequency band 
(41.5-169.25) Hz by the LIGO Scientific Collaboration and 
Virg o Collaboration using data from th e S5 two-year science 
run (|LSC VIRGO Corlaboratio^H 2009) . This limit improved 
on previous indirect limits from the big bang nucleosynthe- 
sis and cosmic microwave background at around 100 Hz. As 
ET will be able to detect GW background signals around 
six orders of magnitude below this limit, it is possible that 
astrophysical GW background signals could form an addi- 
tional 'noise' component concealing the background signal 
from primordial processes. 

The organization of the paper is as follows: In Section 
2 we discuss the rate estimates of SL-GRBs and in Sec- 
tions 3 and 4, we describe in more detail the two previously 
mentioned GW emission mechanisms that could result from 
newly formed magnetars in SL-GRBs. In Section 5 we de- 
scribe how we will calculate a GW background spectrum 
and in Section 6 we discuss issues relevant to detection. In 
Sections 7 and 8 we present our estimations of the GW back- 



http : / /heasarc . gsfc . nasa.gov / docs /swift /swiftsc. html 



LIGO — http: / /www. ligo.caltech.edu/ 
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4 VIRGO — http://www.virgo.infn 

5 ET — http://www.et-gw.eu/ 
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Rate estimate in 




Gpc _3 yr _1 




1800 


r P 


200 




40 



Table 2. The sub-luminous GRB rate estimates used in this 
study. The estimates are denoted by upper, rjj, plausible, rp, 
and lower, n,- 



ground signal from our two single-source emission mecha- 
nisms and finally draw our conclusions in Section 9. 



2 RATE ESTIMATES OF SUB-LUMINOUS 
GRBS 

Table[2]shows rate estimates of SL-GRBs from studies span- 
ning the past four years. Estimates are generally determined 
by statistical arguments, or fits to the log N - log P, peak 
flux, or 'brightness distribution' of bursts. Statistical argu- 
ments are typically based on the two closest sub-luminous 
bursts: GRB 980425 and GRB 060218, detected within 2 
years of operation by Swift. As rates based on Poisson statis- 
tics could be affected by small number statistics, some au- 
thors choose to fit to the log N - log P distribution of ob- 
served bursts. Using this method, a SL-GRB population can 
be accounted for by decreasing the lower bound of the lumi- 
nosity function (LF) or by employing a two-component LF 
i|Cowardll2005[ ). 

The table shows that estimates extend over a range 
(40 — 1800) Gpc _3 yr _1 , reflecting the present uncertainties 
on the nature of these bursts. For example, other than un- 
certainties in the LF, it is still not clear if these bursts are 
LGRBs viewed off-axially o r are an intri nsically different 
population (see discussion in lCowa rd 2005). Assuming that 
SL-GRBs are a unique population with an intrinsic differ- 
ence in central engine from LGRBs, the rates we will adopt 
in this study are shown in Table [2] 

For a plaus i ble ra te we take the most recent estimate 
of IVirgilii et~aD (|200Sh . r P = 200 Gpc~ 3 yr _1 . As shown in 
Table [2] this estimate is of a similar order to the most 
likely values published in the other studies. As an upper 
limit we take the larg est estimate shown in Tab le [5] of 
ru = 1800 Gpc" 3 yr- 1 (|Guetta fe Delia Valid 120071) - this 
value is ~ 9% of the local rate of SNe Ib/c. We note that 
this is a similar fraction of SNe Ib/c producing magnetars 
to that of Type II SNe - around 10% - as sug gested by 
iMurase et ail (|2006T l and ISoderberg et all (|2006h . For our 
lower bound we take a value of tl = 40 Gpc _3 yr _1 . We 
obtain this value by taking a typical lower bound of around 
100 Gpc _3 yr _1 based on the estimates shown in Tableland 
in correspondence with LGRBs observed during the Swift 
Era, we a ssume that 40% o f SL-GRBs will also have X-ray 
plateaus (|Evans et alj|2009h . 



3 TRIAXIAL GW EMISSION FROM 
MAGNETARS 

Rotating NSs with a triaxial shape have a time varying 
quadrupole moment and hence radiate GWs at a frequency, 
/, which is twice the rotational frequency. A NS born with a 
rotational period Po looses rotational energy through mag- 
netic dipole torques and GW emission: 

di?rot _ di?dip . dE e 



+ 



dt dt dt 

with rotational, dipole and gravitational energy loss rates 
d.E ro t 



(1) 



2 T Af 

dt :7rJzz/ d?' 



d£ dip v , 4 tt 4 R 6 B 2 j4 
= h > :>J = — J 



dt 



dE e 



6c 3 

32^ 6 G/ 2 zP 2 /6 



(2) 



dt 8WJ 5c5 

These result in a change in the frequency: 



df 



K„- 



_ ^dip ,3 , 

dt 7T 2 / ZZ 3 + 7T 2 / Z; 



■f 



(3) 



From the equations for d_E gw /d/ and d//d£, we can write 
the emitted GW spectral energy density as follows: 

K 



d|gw =K f {1 + 
where 



7T 2 7 Z 



-f 2 )- 1 with fe [0,2/Po] 



K 



192tt 4 GJ 3 z p 2 
5c 2 i? 6 B 2 



(4) 



(•») 



In this expression R is the radius of the star, p = (/ xx — 
7 yy ) /7 ZZ is the ellipticity in terms of the principal moments of 
inertia, B = B s sin a where B s is the surface polar magnetic 
field strength and a the angle between the rotational and 
dipole axes. 

The majority of NSs are understood to be born 
with magnetic fields of the order of 10 12 — 10 13 G 
and rotational peri ods of the order of tens or hun- 
dreds of milliseconds ([Regimbau fe de Freitas Pacheco|[2000l : 



iFaucher-Giguere fe KaspJ 120061 : ISoria et al.l l2008h and will 
make a negligible contribution. However, NSs with suf- 
ficiently high initial rotational periods, ~ 1-3 ms, which 
undergo violent convection or differential rotation dur- 
ing the first seconds after birth can obtain super-strong 
crustal magnetic fields (B a in the range ~ 10 14 - 10 16 G 



2 



through an q-Q dynamo act ion (|Duncan fe ThompsonlT99 
iThompson fe Duncan! 1 19931 ) . For these highly magnetized 
NSs, the distortion induced by the magnetic torqu e can 
become significant (|Palomball2000l : Ik 011110 et all [2000). and 
GW emission can be orders of magnitudes larger than for 
ordinary NSs l|Palomball200ll ). Because it could carry away 
most of the initial rotational energy of millisecond magne- 
tars, this scenario provides a natural explanation for the 
absence of the signature of enhanced energy injections in X- 
ray spectra of superno va remnants around known magnetars 
l|DallOsso et al.ll2009h . 

In this study, we take Po = 1 ms, and corresponding to a 
typical NS of mass 1AM P) take R = 10 km and J zz = 1 x 10 3S 
kg m 2 l|Arnett fe Bowers! Il976l : iBonazzola fe Gourgoulhonl 
1996). We consider two different scenarios corresponding to 
different configurations of the magnetic field: 
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Table 1. A sample of the published estimates on the rate of SL-GRBs along with a brief description of how 
the rate was determined. 



Reference 


Rate estimate 
Gpc _3 yr _1 


Notes 


Soderbere et al. ( 20061 


230±f™ 


A Poisson statistical estimate based on the detection volumes for 
bursts similar to GRB 980425 and XRF 060218. 


Pian et al. C20061 


110 +180 


A fit to the log TV - log P distribution of BATSEj^ data using 

a smoothed broken power law LF with a lower bound set by GRB 980425. 


Guetta & Delia Valle (20071 


3801^06 

200-1800^ 

110-1200£ 


Poisson statistical estimate determined as well as two fits to the 

log N — log P distributions of both Swift and BATSE. A single power law 

LF was used with a lower bound based on GRB 980425. 


Liang et al. (20071 




The LF and rate density are estimated using Swift bursts with known z. 


Chapman et al. (20071 


700+ 360 
' uu _360 


Estimates obtained by correlating galaxies within 155 Mpc to BATSE 
bursts with properties similar to known SL-GRBs. 


Virsrilii et al. f 20081 


200+ 200 


LF parameters and 2 values estimated through simulation. Rate estimates 
obtained through statistical comparison with the observed Swift 
luminosity-z distribution. 



a BATSE - the Burst and Transient Source Experiment on the Compton Gamma-Ray Observatory, launched 
in 1991, recorded 2704 GRBs during its 9 years of operation. 
^ Estimate based on Poisson statistics. 
c Estimate based on BATSE data. 
d Estimate based on Swift data. 



(i) Poloidal field For the case in which the inter- 
nal magnetic field is purely poloidal and matches the 
dipolar field in the exter i or, t h e ellipticity is give n by 
iBonazzola fc Gourgoulhonl (|l996l ); iKonno et all (|2000l ): 



Pb =9 



R B 



(6) 



According to the n umeri cal simulations of 
IBonazzola fe Gourgoulhonl (|l996j ). the magnetic dis- 
tortion parameter g of a typical NS, which depends on 
both the equation of state and the magnetic field geometry, 
can range from 1 — 10 for a non-superconducting interior 
and can increase to 100 — 1000 for a type I superconductor 
in which all the magnetic field has been expelled from 
the superconducting core. It can take on an even greater 
values for type II supercon ducting cores or counter rotating 
electri c currents. Following iRegimbau fc de Freitas Pacheccl 
(2006) we take g = 520 in the mid-range of permissible 
values of a type I superconducting core, corresponding to a 
core whose dimension represents 95% of the equator radius. 
We take B = 5 x 10 14 G as representative of the magnetar 
population, based on average observational values of soft 
gamma repeaters and anomalous X-ray pulsarijf], but note 
that we have excluded two objects with characteristic times 
of 230 kyr and > 1300 kyr for which dissipation of the 
magnetic field may have been important. 

We find ellipticity ps = 4.8 x 10 -4 . To have an idea of the 
magnitude of such deformation, we can compare this elliptic- 
ity with the ones predicted for elastic deformations of NSs. 



6 The McGill SGR/AXP online catalog can be found at 
http://www.physics.mcgill.ca/~pulsar/magnetar/main.html 



In this respect, ps = 4.8 x 10~ 4 is about two orders of magni- 
tude larger than the ma ximum elastic q uadrupole deforma- 
tion of conventional NSs ((Horowitz 2010), but comparable to 
elastic deformations sustainable by solid strange stars, and 
1 — 2 orders of magnitude below the upper limi t derived for 
crystallin e color-superconducting quark matter |Owenll2005l ; 
lLhill2007l ). For p B = 4.8 x 10~ 4 , the GW emission is negligi- 
ble compared to the magnetic torque and eq. ([4]) simplifies 
to (K « n 2 Izzf- 2 ): 



dE s 



df 



Kf 



We note that for a rotation period of the order of ms, 
strongly magnetized relativistic winds could slow down the 
star in a f ew minutes as energy i s rapidly transferred to 
the ejecta ([Thompson et al. 1 12004 iBucciantini" et al. 1 120071 : 
iMetzger et al.l 2007) . However, this effect is expected to be 
negligible for B < (6 - 7) x 10 14 G. 

In theory, one could consider values up to g = 1000 — 10000 
and B = 10 16 G for which GW emission would dominate 
in most of our frequency range (K >> it 2 I zz f~ 2 ). How- 
ever, this scenario would produce ellipticities of order unity, 
much higher than the ones typical ly considere d for magnetic 
|DallOsso et al.ll2009l) or elastic |Owenll2005h NS deforma- 
tions. 

(ii) Toroidal field A number of studies have suggested 
that the internal magnetic field B t could be dominated 
by a strong toroidal component in the range 10 15 — 10 17 
G llCutlerll200l: Stell a et al ] |2005l : IPallOsso fc Stellall2007i : 
IPallOsso et al.ll2009h . which could induce a prolate defor- 
mation with ellipticity 
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p B ~ 1.6 x 10" 



B t 
10 16 



(8) 



Following IPallOsso et all (|2009l ). we assume a pure internal 
toroidal field of B t = 10 16 G , with an extern al magnetic 
field of the order 10 14 G. B t was deduced by IStella et al.l 
<|2005l ) from studies of the energetics and likely recurrence 
time of the 2004 December 27 event from SGR 1806-20, 
and is consistent with the thermal emissions observed in 
Anomalous X-ray Pulsars, ass uming they are power ed by 
the decay of the magnetic field jKaminder et al.ll2007l ) . Thi s 
value also supports constraints set bv lVink fe Kuiperl (|2006T ) 
on the X-ray spectra of the supernova remnants surrounding 
known magnet ars. 

In this case, both GW and magnetic dipole losses con- 
tribute to the magnetar spindown. At small frequencies 
(/ J; 100 Hz), the emission is dominated by the electro- 
magnetic contribution, but above / ~ 500 Hz, GW emis- 
sion becomes the most important process, approaching the 
saturation regime where spindown is purely gravitational 
{K » 7T 2 / ZZ /- 2 ): 



dE s 



df 



(9) 



Increasing B t to 5 x 10 16 G or 10 17 G lowers the frequency at 
which GW emission becomes the dominant contribution to 
around 100 Hz or 30 Hz. 



4 GW EMISSION FROM SECULAR BAR 
MODE INSTABILITIES 

Bar-mode instabilities associated with NS formation derive 
their name from the 'bar-like' deformation they induce, 
transforming a spheroidal body into an elongated bar 
that tumbles end-over-end. The highly non-axisymmetric 
structure resulting from a compact astrophysical object un- 
dergoing this instability makes such an obj ect a potentially 
stron g source of grav it ation al radiation JChandrasekhaj 
19691; IChandrasekharj Il97d; iFried man fc Schuta 11978 



Lai fc Shapiro! " 11995 
Shibata et aiT 200C 



Brownl |2000|; [New et all 1200 



Ou et al.l 



2004 



Shibata fc Karin 



2004 ; iBaiotti et all [20071 ; iDimmelmeier et all 120081 ; lo"tt 
20091 ). 

A system susceptible to bar-mode deformation is 
parametrized by the stability parameter, P — T/\W\, where 
T is the rotational kinetic energy and W is the gravitational 
potential energy. 

There exist two different timescales and mechanisms for 
these instabilities. Uniformly rotating, incompressible stars 
are secularly unstable if P > 0.14, and have a growth time 
that is determined by the time-scale of dissipative processes 
in the system (such as viscosity or gravitational radiation) 
- usually much longer than the dynamical time-sc ale of the 
system (see e.g. ISaiio et al.ll200ll ; [New et al.ll2000l V In con- 
trast, a dynamical instability sets in when /3 > 0.27, and 
has a growth ti me of the order of the rotation period of the 
object (see e.g. iNew et al.|[2000t ). This is expected to be the 
fastest growing mode. 

Because bar-mode instabilities are a potentially im- 
portant source of gravitational radiati on, they have been 
the subject of many numerical studies. IDimmelmeier et al.l 




100 120 
frequency (Hz) 



Figure 1. The rest frame energy spectrum for GW emission by 
a secular bar mode instability for parameter values M = 1.4Mq; 
R = 20 km; B = 10 14 G; P = 0.20. During the first 2000 s of the 
evolution, the GW energy is emitted in the 50-150 Hz energy 
band. 



(2008) found that the post-bounce core cannot reach suf- 
ficiently rapid rotation to become unstable to the classical 
high-/3 dynamical bar-mode instability. However, they found 
that many of their post-bounce core models had sufficiently 
rapid rotation to beco me subject to a low-/3 d ynamical insta- 
bility first by seen bv lCentrella et all (|200ll ). The potential 
for enhancements in the GW emissions by dynamical insta- 
bilities at low /3 is encouragi ng and has been demonstrated 
in a number of other studies (jshibata et a.1 .1120021 : lOtt et aD 
120051 ; IScheidegger et al 1 120081 . |2010| V 

The requirement of rapid rotation for post-collapse in- 
stabilities suggest that a GR B progenitor, typ ically required 
to be in high rotation ([Wooslev fc Jan ka 2005), may provide 
favourable condi tions for such instabilities to set in. In this 
paper we follow ICorsi fc Meszarosl (|2009al lbl) and consider 
GW emissions from the longer lived secular bar-mode insta- 
bility possibly associated with the observed shallow decay 
phase observed in GRB X-ray afterglows discussed earlier. 
In the next section we describe the single source models we 
employ to estimate the GW backgrounds from this instabil- 
ity. 



4.1 Single-source spectrum from secular 
bar-mode instabilities 

The single-source emission mechanism used in this se ction 
is motivated by the study of lCorsi fc Me szaros (2009a) who 
estimated the GW emissions from a secular instability in 
a newly born magne tar. Their work extended the work of 
lLai fc Shapiro! (| 19951 ) for the quasi-static evolution of NSs 
under gravitational radiation. Treating the NS as a poly- 
trope of index n = 1, they assumed typical parameter values: 
total mass, M = 1.4M©; radius, R = 20 km; initial magnetic 
dipole field strength at the poles, B = 10 14 G; and /3 = 0.20 
corresponding with the middle of the expected range for the 
secular instability (0.14 < P < 0.27). They estimated quasi- 
periodic GW emissions at around 150 Hz, with characteristic 
amplitudes h c ~ 10 -21 at 10 Mpc. 

Figure [1] shows the rest frame energy spectrum, 
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d-Bg W /d/, from this mechanism. This function was computed 
through 



dE s 



df 



dt 



dt 
df 



(10) 



using data for the lumino sity and frequency evolut ion of 
the instability supplied by ICorsi fc Meszaro s (2009a). The 
bulk of the emission takes place over the first 2000 s; dur- 
ing the first 100 s the GW frequency is constant at 150 Hz. 
This is then followed by a slow decline over about 2000 s to 
~ 100 Hz. Their analytical treatment included the effects 
of energy losses from magnetic dipole radiation. In princi- 
ple, the secular evolution would bring the star to reach a 
stationary football configuration; thus one could follow the 
predicted GW signal unt il its frequency approaches zero. 
ICorsi fc Meszaros! (|2009aT ) conservatively shut off the bar 
emission after about few thousand seconds of evolution (as 
for typical GRB plateau durations), when the GW signal 
had a frequency of about 50 Hz, and its amplitude was falling 
below the ALIGO sensitivity curve. We note that the ET 
sensitivity below 50 Hz is one order of magnitude better than 
for ALIGO, but the final stages of the secular evolution are 
indeed highly uncertain (due to e.g. viscosity effects possi- 
bly coming into play, see e.g. Lai and Shapiro 1995, Corsi 
and Meszaros 2009) . Thus, also in this analysis we conserva- 
tively assume that the bar emission completely shuts off on 
the typical duration of GRB plateaus, when the GW signal 
frequency is around 50 Hz. 

We note here that, despite uncertainties on the values 
of the model parameters, FigfTJcan be considered as an av- 
erage case: mass, radius and magnetic field of the NS are 
chosen so as to represent the typical case for a newborn 
magnetar; /3, is chosen to be in the middle of the secular 
instability range. We also stress that, as discussed in Corsi 
& Meszaros 2009, with this typical choice of parameters, the 
observed timescale of GRB plateaus is correctly reproduced. 
Moreover, the amount of energy released by the magnetar 
for dipole losses during such timescale, is of the order of 10 50 
ergs, i.e. comparable to the isotropic energy output of long 
sub-luminous GRBs, and thus sufficient to actually cause a 
visible plateau in their light curve. 



5 THE GRAVITATIONAL WAVE 
BACKGROUND 

5.1 Cosmic metallicity evolution 

The spectral form of an astrophysical GW background 
(AGB) is highly dependent on the variation of the 
event rate with z jHowell et all 12004 ICoward et"afl l200ll ; 
iRegimbau fc Mandidl200Sl ). In general, due to the relative 
short lifetimes of massive stars (of order tens of Myr) the 
transient populations that produce AGBs are assumed to 
track the star formation rate (SFR); for transient popula- 
tions of coalescing compact objects an additional factor is 
included to account for merger time delay. For the case of 
an AGB from SL-GRBs, there is growing evidence th at cos- 
mic metallicity evo lution must also be considered (|Lill2008l ; 
iModiaz et alfeosft . 

For WR stars to retain sufficient rotation to power 
a GRB, angular momentum losses thr ough stellar-wind 
induced mass-loss must be minimized |Wooslev fc Hegerl 



2006). As wind-driven mass loss in WR stars is under- 
stood to be dependent on a high enough fraction of iron, 
a low-metallicity environment is an essential requirement 
jVink fc de Koteij|2005l ; IWooslev fc Jankall2005» . 

To account for metallicity evolutio n of SL-GRBs with 
redshi ft we adopt the simple model of ILanger fc NormanI 
l|2006h 



r(0.84,e 2 10°' 3z ) 
r(0.84) 



(11) 



where e = Z/Zq is the fraction of solar metallicity and 
&(z, e) is the fraction of massive stars at z born with metal- 
licity less than Zq6. Here, T a nd V are the incomplete 
and complete gamma functions. ILanger fc NormanI (2006) 
found that a metallicity cutoff of e = 0.1, corresponding 
to ~ 1 GRB per 100 WR stars throughout the Universe, 
was able to reproduce the observed global ratio of the rates 
of GRBs to core-collaps e SNe . This value was also used by 
ISalvaterra fc Chincarinil (2007) who suggested that luminos- 
ity evolution was required to reproduce the Swift distribu- 
tio n at high z. A n analy sis of the Swift data (to August 2009) 
by iButler et al. (2010) ruled out luminosity evolution and 
found a more relaxed cutoff e ~ 0.2 — 0.5, was adequate to 
reproduce the observed sample. Their more modest depen- 
dence on metallicity was supported by studies of the mas s 
distribution of GRB host galaxies bv lKocevski et aD {2009). 
Based on these studies, we adopt here a range e ~ 0.1 — 0.5 
to allow for present uncertainty. 



5.2 Source rate density evolution 

Our source rate evolution model for SL-GRBs with redshift, 
Rsl(z), is obtained by scaling the star formation historjQ, 
Rsf(z), with the function ty(z,e): 

R SL (z) = ^(z,e)R S F(z). (12) 

For R sf(z), we employ the model of IHopkins fc BeacomI 
( 2006), who constrained the star formation rate history by 
combining recent measurements taken from sources observed 
at ultraviolet, far-infrared and radio with previous more ro- 
bust measurements taken over the last decade. By normal- 
izing Rsl(z) to the the local (z — 0) rate, we produce a 
dimensionless evolution factor 



e(z) = R S L(z)/RsL(z = Q). 



(13) 



This allows us to extrapolate a local rate density to cosmo- 
logical volumes. 



5.3 The event rate equation 

In order to evaluate the contribution by a population of GW 
sources to a stochastic background signal, knowledge of the 
rate of emissions from recent and past epochs is essential. 
The source rate evolution of an AGB is modeled by the 
differential event rate, dR/dz, which describes the rate of 
events in the redshift shell z to 2 + dz: 

dV r e(z) 
dR = — : — — dz . 



dz 1 + z 



(14) 



7 In units of mass converted to stars per unit time and volume. 
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Figure 2. The differential event rate dR/dz of SL-GRBs 
as a function of redsh ift based on the SFR model of 
iHopk ills fc Beacoml J2006h . The three curves correspond to dif- 
ferent metallicity cutoffs e = (0.1; 0.5; ^)Zq at the upper rate 
estimate of ry = 1800 Gpc -3 yr -1 . Inclusion of cosmic metallic- 
ity evolution increases the contribution from sources at higher z. 
This is shown by a shift in the peak of dR/dz from z ~ 2 (e = 1) 
to z ~ 3 assuming SL-GRBs follow a low metallicity dependence 
of e = 0.1Z e . 



The (1 + z) factor accounts for the time dilation of the 
observed rate by cosmic expansion; its inclusion converts 
source-count information to an event rate. The parameter 
ro is the local rate density, usually defined within a vol- 
ume spanning the Virgo cluster of galaxies (at around 20 
Mpc). This factor is fundamental to estimating the number 
of potentially observable events and is determined using es- 
timated source rates within a larger fixed volume of space. 
For the factor ro, we adopt the values discussed in section 2 
of (ru,r P ,r L ) = (1800,200,40) Cpc'V" 1 - 

The co-moving volume element dV describes how the 
number densities of non-evolving objects locked into Hub- 
ble flow are constant with redshift. This is obtained by c al- 
culating the luminosity distance from (cf. Peebles 199^, p. 
332) 



■ , n f i n c f z dz' 

d L (z) = i + *hr / T7 jt j 

H o Jo h \ z ) 
and using eq. 3 ofH orciam fc Madaul (|2000h . 

dV _ 4irc d L 2 (z) 
dz ~ H (l + z) 2 h(z) ' 

The normalized Hubble parameter, h(z), is given by 

h{z) = H{z)/H = [n„(l + zf + Q A ] V2 , 

for a 'flat-A' cosmology. We take Sl m = 0.3, = 0.7 and 
Ho = 70 km s - 1 Mpc -1 for th e Hubble parameter at the 
present epoch (Rao et al, 200(j). 

Figure [5] shows the form of dR/dz assuming different 
metallicity cutoffs for SF-GRBs. In this plot we compare 
curves corresponding to the different cutoffs e = (0.1,0.5) 
with the e = 1 solar metallicity curve. The inclusion of 
cosmic metallicity evolution increases the contribution from 
higher-z sources, shown by a shift of the peak from z ~ 2 



(15) 



(16) 



(17) 



to z ~ 3 for a low metallicity dependence, e = O.IZq. The 
magnitude of the peak value also increases by up to a fac- 
tor of 3. Thus, the additional contribution from sources at 
higher z should enhance the lower frequency component of 
the stochastic background signal through redshift. 

5.4 The Gravitational Wave background spectra 

The spectral time- integrated flux density or spectral fluence, 
in J m -2 Hz -1 , of a quadrupole GW signal at a luminosity 
distance dh(z) from a single source can be expressed as 



Fss ( /obs , Z) — 



d£ gw (/ obs ) (1 + 2) 



(18) 



df 47Td L (z) 2 ' 

where d-E gw (/obs) is the spectral GW energy at the observed 
frequency / b s , which is related to the source frequency / by 
the redshift factor: / b s = //(l + z). 

The background spectral flux density, in W m~ 2 Hz -1 , 
from all events throughout the Universe is obtained by in- 
tegrating the product F BB (f ohB , z)dR/dz over the redshift 
range z = to zn m : 

Fs(/obs) = / hm [F ss (f ohs ,z)(dR/dz)] dz, (19) 
Jo 

with the integrand given by (|14[) and (|18[1 within a limiting 
redshift, z\i m = 10, which we take as the beginning of stellar 
activity. In support of this value we note that GRB 090423, 
the most distant recorded burst (z = 8.2), showed no evi- 
dence of properties that were inco nsistent with the m ajority 
of the observed GRB population (|Tanvir et a.lj[200sh . 

The spectral energy density of a GW background is con- 
ventionally expressed by the dimensionless energy density 
parameter, fiB(/obs), defined as the energy density of GWs 
per logarithmic frequency interval normalized to the cosmo- 
logical critical energy density p c = 3Hq / 8irG. This func- 
tion can be constructed from Fsif^) (|Ferrari et aTH l999a; 
iHowell et al.ll2004 iRegimbau fc Mandidl2008l )~ 



fWobs) =/F S (/ obs )/(p c C 3 



(20) 



Throughout this paper we will present our estimated GW 
background spectra using this function. 

5.5 The duty cycle of an astrophysical GW 
background 

When we consider AGB signals, in addition to the energy 
density and characteristic frequency, another useful quan- 
tity is the duty cycle. DC dBlair fc Jull 19961; lMaggiorejl2000l ; 
ICoward fc Regimbaull2006l ; IRegimbau fc Mandicfeoosi ) . The 
value of the DC is given by the ratio of the typical dura- 
tion of the signal, r, to the average waiting time between 
the reception of successive events. The waiting time will de- 
pend on the rate of events as observed in our frame, R, and 
thus a DC is generally defined by the quantity Rxt. When 
considering a cosmological distribution of events, the DC is 
determined by sources within a limiting redshift, zu m : 



DC{z lim ) 



(1 + z)r (dR/dz) dz , 



(21) 



Here, the signal duration is dilated to (1 + z)t by the cos- 
mic expansion and the quantity di?/dz is the cosmologically 
dependent differential event rate given by equation 1141 
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Many studies are concerned with the value of DC as 
seen at the detector and determine the value of DC by set- 
ting zii m equal to the redshift at which stellar activity began. 
In this case equation (|21|) provides a total value, DC(z\\ rn ). 
As source rate evolution will increase out to large cosmolog- 
ical volumes, it is interesting to see how DC too increases 
with z. In this study, we will therefore calculate DC as a 
function of redshift. 



In general, for an AGB the signal is defined as 
continuous if it has a DC of unity or above. As it 
will still be possible to resolve individual events in this 
regime, a more conservative threshold can be defined 
from the view of s i ngle detections as DC ^ 10 (see 
iReeimbau fc Hughes! 120091). Thus, using the convention 
of Regimbau fc de Freitas Pachecd (|200fjT ). non-continuous 



signals can be further categorized into shot and popcorn 
type. More descriptive definitions of these components are 
provided as follows: 




10 



10 10 
frequency (Hz) 



Figure 3. The design sensitivity curves for advanced LIGO 
(ALIGO), Advanced Virgo (AVIRGO) and two possible config- 
urations of the third generation Einstein Telescope (ET): ET-B 
and ET-C. See Sec. 7 for more details. 



Continuous (DC ^ 10): This signal is, at any given 
time, the superposition of a large, but random, number of 
overlapping signals. As the amplitude of each contributing 
signal is itself random, the central limit theorem will apply 
leading to a Gaussian distribution in amplitudes. Because 
it will be difficult to resolve the individual components, 
this signal can p otentially mask a background signal of 
primordial origin (jMaggiorel 120001 ). For this reason, AGB 
with DC ^ 10 would be bad news from the perspective of 
primordial background searches. 



Popcorn noise (0.1 ^ DC < 10): This signal will 
manifest in GW data as a non-continuous stochastic back- 
ground signal with an amplitude distribution dependent 
on the spatial distribution of the sources. For background 
signals with DCs at the lower end of the popcorn noise 
range, the individual components will rarely overlap. In this 
case the signal will be dominated by the closest events of a 
population of sources. 



Shot noise DC < 0.1 : For this signal, the waiting 
time between ev ents is large in comparison with dura tion 
of a single event (|Regimbau fc de Freitas Pachecd [200^ 1. 



Using these thre e definitions of DC, 
ICoward fc Regimbau! (2006) have shown that an AGB 
can be divided into three different detection regimes, each 
defined by the corresponding shells of zn m . For most types 
of AGB, at low z, the signal produced is predominantly of 
the shot-noise type, extending to popcorn and continuous 
with z, due to time dilation and increased source rate 
density. The weighting of the different AGB regimes for a 
particular source population depends strongly on the event 
rate and is an important consideration when selecting the 
most appropriate signal detection strategy. 



6 DETECTION 

In this study we consider the design sensitivities of second 
generation interferometric detectors, such as ALIGC0 (ex- 
pected online in 2015) and Virgqjl 10 !. and third generation 
ones, for which we will use ETrnfpossible construction will 
be late in the next decade). 

For ALIGO, we use the sensitivity curve based on the 
zero detuning, high laser power configuratiorT 2 !. For ET 
we consider two target sensitivities. Firstly, ET-B, which 
is an underground based design, incorporating lo ng suspen- 
sions , cryogenics and signal and power recycling (Hil d et al.l 
2008). Secondly, a so-called Xylophone configuration, ET- 
C, which merges the output of two detectors specializing in 
different frequency bands: a) an underground low- frequency 
cryogenic configuration with long suspensions and moder- 
ate laser power; b) a high frequency detector implementing 
squee zed light states, large test masses and a high power 
laser (|Hild et al.|[2010l ). The advantage of this strategy is 
that it can decouple the obstacles in operating a high power 
laser alongside a c ryogenically coole d suspension optimized 
for thermal noise (|Shoemakerll200ll l. The design sensitivity 
curves for these four detectors are shown in Figure [3] 

The most promising detection strategy for continu- 
ous GW background signals is cross-correlating the out- 

8 http://iopscience.iop.Org/0264-9381/27/8/084006 

9 https:/ /pub3.ego-gw.it/codifier/includes/showTmpFile.php?doc 
=2219&calledFile= VIR-027A-09.pdf 

10 www.virgo.infn.it 

11 A three year design study for the Einstein 
GW telescope began in May 2008. For details see 
|http://www.ego-gw.it/ILIAS-GW/FP7-DS/fp7-DS.htm| De- 
sign Study Proposal for E.T. (Einstein Gravitational Wave 
Telescope), submitted to the EU Seventh Framework Pro- 
gramme. 

12 The ALIGO sensitivity curve is described 
in the public LIGO document |ligo-t09002"88l v2 
(https://dcc.ligo.org/public/0002/T0900288/002/AdvLIGO 
noise curves.pdf). 
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put of two neighbour ing detectors (see, iMaggiorel l200d : 
I Allen fc Romano) 1 19991) . For this strategy to be achieved, 
the detectors must be separated by less than one reduced 
wavelength, which is about 100 km for frequencies around 
500 Hz where Qsif) might peak. The detectors also need 
to be sufficiently well separated that their noise sources are 
largely uncorrelated. We note that although this may not be 
possible for ET, techniques are in development to remove en- 
viron mental noise and instrumental correlations l)Fotopoulosl 
120081 ). 

Under these conditions, assuming Gaussian noise in 
each detector and optimal filtering, a filter function cho- 
sen to maximize the s ignal-to-noise ratio, SN R, for two such 
detectors is given by (|Allen fe Romano! 1 19991 . eq.3.75) 



50tt 4 



T 



7 2 (/)"l(/) 
f a S nl (f)S n2 (f) 



df- 



(22) 
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Here 7(/) is the 'overlap reduction function', which accounts 
for the separation and relative orientation of the detectors, 
and Sni(f) and 5^2 (/) are the noise power spectral densities 
of the detectors; T is the integration time. As the optimal 
filter depends on a range of filter functions based 

on theoretical expectations of this function will need to be 
used. In this study we adopt a value of SNR — 3 to indicate 
detection, corresponding w ith false alarm rate of 10% and 
detection rate of 90% and (| Allen fc Romano! [l999h . 

For signals of the shot noise and popcorn type, the 
standard cross correlation strategy in the frequency domain 
may not be optimal and other methods in the time do- 
main have been proposed or are curre ntly under investiga- 
tion i n the LIGO- Virgo collab oration (|Drasco fc Flanagan! 
120031 ; ICoward fc Burmarj|2005l ) . For signals of the shot noise 
type, individual events may be clearly distinguishable, and if 
within a detector's range, can potentially be resolved using 
data analysis techniques such as matched filtering. 

To estimate the detectability of the GW backgrounds 
considered in this paper we will assume continuous signals, 
and hence use the cross correlation statistic to determine 
SNRs. For cases in which a significant proportion of the 
background signal will not be continuous, we will also in- 
vestigate the different z regimes in which the shot, popcorn 
and continuous regimes exist. We calculate the SNRs for the 
two generations of GW detector outlined above; for second 
generation we assume 3 years of integration by an ALIGO 
configuration; for third generation we assume 1 year of inte- 
gration by (ET-B; ET-C). We further assume: a) separated 
detectors; b) an optimal case in which a pair of equivalent 
detectors are situated within several km and aligned. For 
ALIGO we will employ the LIGO Hanford/Livingston pair 
(Hl-Ll ) for scenari o a), u sing for 7(/) the form given by eq. 
3.26 of lAllen et all |2002). For ET we assume two detectors 
located in Cascina, of triangular shape (60° between the two 
arms) and separated by an angle of 120°. In the frequency 
range we are interested in (1-1000 Hz), 7(/) reduces to a 
value of -3/8. For case b) we will assume 7(/) = 1 for both 
second and third generation detectors. 

For convenient comparison of £2b(/) with the GW in- 
terferometric sensitivity curves discussed above, the noise 
power spectrum, S n (f) (in units of Hz -1 ) over a frequency 
range A/ can be expressed in terms of a detector energy 
density, Qdet(f), over an integration time Tint: 



Figure 4. The GW background spectrum from triaxial de- 
formations in newly born magnetars associated with SL-GRBs 
in the regime where the spindown is dominated by the mag- 
netic torque. The three curves assume that triaxial deforma- 
tions are introduced by an internal poloidal magnetic field 
(B = 5 X 10 14 G;g = 520; / zz = 10 45 kgm 2 ;R=10km). The 
curves are presented for three rates of occurrence: (»"u j r P j r L ) = 
(1800, 200, 40) Gpc _3 yr~ 1 and metallicity cutoffs of e = 0.1 (thick 
curves) and. e = 0.5 (thin curves). The sensitivity curves of sec- 
ond and third generation laser interferometric detectors are rep- 
resented by ALIGO and ET-B, in terms of fidet(/) assuming a 
3 year integration. Based on observational values of SGRs and 
anomolous X-ray pulsars, we assume the value of B used here 
is representative of the magnetar population. This value gives 
pg = 4.8 X 10 — 4 , comparable to elastic deformations sustainable 
by solid strange stars, and 1-2 orders of magnitude below the 
upper limit derived for crystalline color-superconducting quark 
matter. 



fidet(/) 



507T 2 



Sn(f) 



ml 



A/r int7 (/) 



(23) 



The sensitivity curves will be presented for the optimal sce- 
nario b), as discussed above for integration times of 1 and 3 
years for ET and ALIGO respectively. 



7 THE GW BACKGROUND FROM NEWLY 
FORMED MAGNETARS 

Figures|3]to|S]show the function fiB(/obs) from triaxial defor- 
mations in magnetars associated with SL-GRBs for the three 
mechanisms discussed in section [3] Curves are displayed at 
the rates (ru,rp,?~L) and metallicity cutoffs e = 0.1 (thick 
lines) and e = 0.5 (thin lines). Each plot also includes the 
detector sensitivities of ALIGO and ET-B assuming 3 years 
of integration. 

Figure [4] shows the function r2p,(/obs) assuming that 
the internal magnetic field is purely poloidal and matches 
to the external dipole field. For this case the GW emission 
is negligible in comparison with the magnetic torque and 
the gravitational signal increases as f until a maximum at 
around 840 Hz for rate nj. We see that the signal is outside 
the sensitivity of ET, and even at the most optimistic rates 
this background would not be detected within a reasonable 
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Figure 5. As for Figure[4] but in the regime where both magnetic 
and GW emission contribute to the spindown. The triaxial defor- 
mations are produced by an internal toroidal field (B t = 10 16 G; 
B = 10 14 G). 



integration time. For higher values of B or ellipticity and 
smaller values of 7 ZZ , the amplitude increases until the GW 
emission dominates (fi ~ f 2 ) at large frequencies (see fig. 
5), eventually reaching a saturation regime (fig. 6). 

Figure [5] displays the function fiB(/obs) assuming that 
internal toroidal fields contribute towards a prolate distor- 
tion. For this case both GW and magnetic dipole emissions 
contribute to the spindown. We see that for the upper rate 
rjj and e = 0.1 the signal peaks at fiB(/obs) ~5x 1CP 10 at 
around 830 Hz. Even at this upper rate only a small part of 
this signal is within the sensitivity of ET-B. 

In Figure [5] we show the upper limit in which spindown 
is purely gravitational (from equation [9]) . We see that the 
background signal increasing with f 2 and reaching a maxi- 
mum fiB(/obs) ~4x 10~ 8 at around 660 Hz for rate r\j. In 
this case, fi depends linearly on 7 ZZ and is independent of 
the ellipticity. 

As illustrated by the thin lines for each curve, which 
represent e = 0.5, we see that a more relaxed metallicity 
cutoff results in a smaller contribution below the peak, and 
hence, a lower SNR, as will be shown below. 

Tables 3 and 4 display the SNRs for the three 
AGBs considered in this section for metallicity cutoffs of 
e = 0.1 and 0.5 respectively. Estimates are shown for 
rp = 200Gpc _3 yr _1 . Values for the rates (rv;rt,) = 
(1800; 40) Gpc _3 yr _1 can be obtained through the ratios 
(200/1800; 200/40). For a GW background produced by 
toroidal field induced distortions, as illustrated in Fig[5l 
at the upper rate ru a small part of this signal is within 
the sensitivity of ET-B. This produces a SNR of 1.8 for 
3 years of integration. Optimistically, for an upper limit 
AGB from pure GW spindown, the SNRs confirm that de- 
tection would require a third generation detector. Assuming 
a metallicity cutoff of e = 0.1, we find that source rates 
of (133; 349) Gpc~ 3 yr _1 would result in detection by (ET- 
B, ET-C) within 3 years at a SNR of 3. A more relaxed 
metallicity cutoff, e = 0.5, increases the corresponding rates 
required for detection to (210; 552) Gpc _3 yr _1 . As shown 



Emission 


ALIGO 


ET-B 


ET-C 


Mechanism 








Poloidal 


i x icr 5 


0.08 


0.03 


field 


(6.1 x icr 4 ) 


(0.07) 


(0.1) 


Toroidal 


2.5 x 10" 5 


0.2 


0.07 


field 


(1.5 x 10" 3 ) 


(0.5) 


(0.2) 


Pure GW 


0.006 


4.5 


1.7 


Spindown 


(0.04) 


(12.0) 


(4.6) 



Table 3. The SNRs obtained through cross-correlation for 
a GW background of triaxially deformed newly born mag- 
netars associated with SL-GRBs for an event rate of rp = 
200 Gpc~ 3 yr -1 . Values for the other rates considered in this 
study, (ru; rx) = (1800; 40), Gpc - 3 yr _1 , can be obtained through 
the ratios (200/1800; 200/40). A metallicity cutoff of e = 0.1 is 
assumed for SL-GRBs. These tabulated values assume 3 years of 
integration by cross-correlation of data from two detectors. We 
also assume the overlap reduction functions described in section 
5 - values obtained by optimally orientated and co-located detec- 
tors are shown in parentheses. The SNRs obtained from the two 
emission scenarios considered in section 3 are shown: distortions 
induced by poloidal fields and toroidal fields; in addition we show 
upper limits that assume spindown is solely gravitational. We 
note that the pure GW upper limit could increase up to 3 times 
the values shown above, as J zz can be 1 -3 times the canonical 
value used in this study jRudermanlll972h . 




frequency (Hz) 



Figure 6. As for Figure [4] but showing upper limits on the GW 
background spectrum from triaxial deformations in newly born 
magnetars associated with SL-GRBs. The three curves assume 
a pure gravitationally induced spindown as given by equation 
with / zz = 10 45 kgm 2 . We note that in this case, f! depends 
linearly on / zz and is independent of the ellipticity 



by Table 4 we see that ET-C, which is optimised for greater 
sensitivity at low frequency, ^ 20 Hz, does not improve on 
the SNR of ET-B. We note that this final scenario can be re- 
garded only as an upper limit, but we consider it as it could 
allow ET to place constraints on this source population. 

We note that in general, independent of the particu- 
lar mechanism driving the deformation, a NS with an el- 
lipticity pb and an external field B, will emit GWs ac- 
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Emission 


ALIGO 


ET-B 


ET-C 


Mechanism 








Poloidal 


4.2 x 10" 6 


0.04 


0.01 


field 


(3.4 x 1(T 4 ) 


(0.1) 


(0.04) 


Toroidal 


1.1 x 1CT 5 


0.1 


0.04 


field 


(8.5 x 1CT 4 ) 


(0.3) 


(0.1) 


Pure GW 


0.004 


2.8 


1.3 


Spindown 


(0.03) 


(7.6) 


(3.4) 



Table 4. As for Table 4 but with a metallicity cutoff of e = 0.5 
for SL-GRBs. 
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Figure 7. The ps — B plane accessible by ET-B for magnetars 
associated with SL-GRBs. The shaded zones, set by the prod- 
uct (pb/W~ 4 ) (10 14 /.B), show the parameter space that can be 
explored for different values of the rate and metallicity e. 



cording to equations (4) and (5). The SNR for the corre- 
sponding GW background signal will therefore depend on 
the combination pb/B. Treating pB and B as independent 
parameters, and assuming that the NS is born in associa- 
tion with SL-GRB, we can make a statement on detectabil- 
ity that is independent of the actual mechanism causing the 
ellipticity. This is done in Table 5 by computing, for each 
value of the rate {rh,rp ,r\j) and of the metallicity cutoff 
e, the minimum (p B /W~ 4 ) (10 14 /B) required to have a de- 
tection with a given detector configuration. The full param- 
eter space is illustrated in Figure [7] for which the prod- 
uct (ps/10 -4 ) (10 14 /B), shown by the diagonal lines in the 
Pb—B plane, divides the plot into detectibility zones (shown 
by the legend). We see that for a rate of tl even pure GW 
emission is out of reach. However, with rp we could access 
extreme values and with ru a large part of parameter space 
is detectable. 

In the next section we will consider the background 
signal from secular instabilities which occur on a shorter 
timescale than the emissions considered in this section, 
~ 1000 s, corresponding with the X-ray plateaus observed 
in some LGRBs. For this signal some analysis of the DC 
will be important. 



Rate 


e = 0.1 


e = 0.5 


(Gpc-V- 1 ) 






40 






200 


15 


pure GW 


1800 


2 


3 



Table 5. The minimal value of the product 
(Ps/ICT 4 ) (B/IO 14 )" 1 required to obtain a SNR of 3 with 
ET-B over 3 years of observation, for the 3 rates considered in 
this paper and for metallicity cutoffs of e = 0.1 and e = 0.5. An 
omitted value implies that a detection is not expected. 

8 THE GW BACKGROUND FROM SECULAR 
BAR-MODE INSTABILITIES 

8.1 The GW background spectrum 

The SNRs estimated for a background signal from secular 
bar mode instabilities in SL-GRBs are shown in Tables [6] 
and [7| for the rates (ry, rp, r^). The estimates indicate that 
ALIGO will require 3 yrs of integration by optimally orien- 
tated and co-located detectors to reach a SNR of 1 for the 
optimistic rate r\j. For ET-B however, this signal can poten- 
tially be detected with a SNR ^ 3 in the more conservative 
hypothesis for SL-GRB rates rr, and detector performances. 

We stress here that these estimates are based on two 
main assumptions: 

(i) that at least 40% of SL-GRBs are associated with mag- 
netar progenitors undergoing a secular-bar mode instability; 

(ii) that the magnetar's parameters are those adopted to 
calculate the single-source spectrum shown in Fig. 1. 

Regarding (i), we note that there is significant uncer- 
tainty in how often SL-GRBs could be associated with a 
secularly unstable magnetar progenitor. This, in turn, im- 
plies a large uncertainty on our rate estimates. To address 
this problem, we have chosen a wide range of values. 

We note again that our lower rate accounts for the 
fraction of LGR Bs showing X-ray plateaus in the Swift Era 
- around 40% (|Evans et al.| [2009). We suggest that is 
a reasonable estimate also in view of the uncertainty (ii) 
underlined above. In fact, the parameter values adopted by 
ICorsi fc Meszarosi <|2009af ) aimed at explaining the typical 
case of a ~ 1000 s duration plateau observed in a LGRB 
with an energy release similar to those of SL-GRBs. It is 
thus more conservative to assume that the spectrum shown 
in Fig. [T] would be realized only in a fraction of SL-GRBs 
similar to the one of LGRBs showing a plateau (~ 40%). 
We note that fraction could be higher in SL-GRBs, since 
dipole energy injection from a magnetar can more easily 
cause visible plateaus on less energetic GRBs. However, it 
is not yet clear whether or not X-ray plateaus are always 
caused by a magnetar - see e.g. iPanaitescul (|2008l ) for 
alternative explanation. In the light of these uncertainties, 
we consider tl as a safer estimate. 

We calculate that rates of (48, 80) Gpc _3 yr- 1 for (ET- 
B, ET-C) are required to achieve a SNR of 3 for 1 year of 
integration by separated detectors. For our conservative rate 
r L , (ET-B, ET-C) will require (1.4, 4) yrs of integration. 

Figure [8] shows the quantity fiB(/obs) for the rates 
(ru, rp, tl), in comparison with the sensitivity curves for sec- 
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Rate 


ALIGO 


ET-B 


ET-C 


Gpc^yr^ 1 








1800 


1 (4) 


113 (300) 


67 (178) 


200 


0.1 (0.4) 


13 (33) 


8 (20) 


40 


0.02 (0.08) 


3(7) 


2 (4) 



Table 6. The SNRs achievable through cross-correlation of 3 yrs 
of data by ALIGO and 1 yr by ET-B or ET-C, for a GW back- 
ground resulting from secularly unstable magnetars associated 
with SL-GRBs. Values are shown for the rates (?"tj > r P i r L ) = 
(1800, 200,40) Gpc-V" 1 and assume a mctallicity cutoff e = 
0.1. The SNR estimates assume the overlap reduction functions 
described in section[6]- values obtainable by optimally orientated 
and co-located detectors are shown in parentheses. 



Rate 


ALIGO 


ET-B 


ET-C 


Gpc~ 3 yr" 1 








1800 


0.7 (3) 


96 (257) 


55 (148) 


200 


0.09 (0.3) 


11 (29) 


6 (16) 


10 


0.02 (0.07) 


2 (6) 


1 (3) 



Table 7. As for Table [6] but with a metallicity cutoff of e = 0.5 
for SL-GRBs. 

ond and third generation GW interferometric detectors rep- 
resented by ALIGO and ET-B. The stochastic background 
signal has a frequency bandwidth 5-150 Hz with a peak of 
^B(/obs) ~ 1CP 9 at around 80 Hz. The thin lines in the 
figure show the function f2B(/obs) assuming a more relaxed 
allowance for metallicity, e = 0.5. This illustrates once again 
how a lower metallicity cutoff results in a greater contribu- 
tion of f2B(/oba) at lower frequency. In Appendix A we will 
further discuss the effect of cosmic metallicity on the GW 
background signal. 

Figure [9] compares the function fiB(/obs) with ET-C. 
In comparison with Fig. [8] we see that at the conservative 
rate, r^, the only significant contribution from this signal is 
from > 50 Hz. As can be seen by a comparison of Tables 
5 and 6, this results in a SNR of around a factor of 2 less. 
At the plausible rate rp there is no contribution below ~ 15 
Hz. Therefore, for the most probable rate estimates, ET-C 
could still pursue a primordial GW background signal at its 
most sensitive frequency bandwidth. 

For ET-B this signal occurs in the most sensitive fre- 
quency regime; a largely continuous signal could therefore 
mask any primordial GW background signal with ^ 
10" 9 within the bandwidth 20-100 Hz. A DC analysis will 
indicate what proportion of sources will contribute to a 
continuous signal. This will have implications for both for 
stochastic searches and for single detections by ET, for which 
the detection horizon may enter the DC ^ 10 confusion lim- 
ited regime. We shall investigate this in the next section. 

8.2 The duty cycle 

For a GW from triaxial deformations in magnetars asso- 
ciated with SL-GRBs we find that for rates (?~tj, rp, = 
(1800, 200, 40) Gpc _3 yr _1 sources outside a volume defined 
by 2 ~ (0.07, 0.1, 0.2) contribute to a continuous signal. This 




frequency (Hz) 



Figure 8. Upper limits on the GW background spectrum from 
secularly unstable magnetars associated with under-luminous 
GRBs for the rates {r v ,r P ,r h ) = (1800, 200, 40) Gpc" 3 yr- 1 . 
The sensitivity curves of (ALIGO, ET-B) are shown in terms 
of 0(j et (/) and assume (3, l)yrs of integration and optimally 
orientated and co-located detectors. The thin curves show the 
background estimates for a more relaxed metallicity dependence, 
e = 0.5. 
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Figure 9. As for Figure|8] but showing an alternative configura- 
tion for the Einstein Telescope, ET-C. 

is the result of a long duration, r ~ 10 6 s in equation 1211 
|Stella et al.ll2005l ). 

Figure [10] shows the duty cycle as a function of red- 
shift for a GW background resulting from secularly un- 
stable magnetars occurring in SL-GRBs. For r, we use a 
value of 1000 s - this approximates to the typical duration 
of an X-ray plateau for a GRB. The plot shows that as 
the rate decreases, the continuous contribution to the back- 
ground is from sources at greater distances. For the rates 
(ru, r P , r L ) = (1800, 200, 40) Cpc'V" 1 we find a DC ^ 10 
is reached at around z ~ (0.5, 1.0, 1.6). Therefore sources 
outside volumes defined by these z values will contribute to 
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1000 




redshift, z 



Figure 10. Duty cycle as a function of redshift for a GW back- 
ground from secularly unstable magnetars in SL-GRBs occuring 
at rates of (40, 200, 1800) Gpc^yr" 1 . The shaded areas show 
three zones of an AGB corresponding to different regimes of DC: 
continuous ( DC ^ 10), popcorn (0,1 ^ D C < 10) and shot noise 
(DC < 0.1) llCoward fe Regimbaul l2006h . We see that sources 
beyond (z = 0.5, 1.0, 1.6) contribute to a continuous signal. 



a continuous background signal. Curves for a more relaxed 
cutoff e — 0.5 are shown by the thinner lines. Referring to 
Figure [2] we see that the effect of metallicity dependence is 
small within z ~ 1. This reflected in the curves of Figure [TU] 

The optimal a nd isotropic (orientation av e raged) hori- 
zon distances (see iReeimbau fc Hughes! 120091 : iDalal et"aD 
120061 . for further definitions) are greatest for the ET-B de- 
tector, at distances of z — (0.2, 0.12), both less than the red- 
shift range in which the signal becomes continuous. Thus, a 
confusion-limited background will not affect the resolution 
of individual sources. 

Figure [11] shows again the results of Fig. [8] but with 
thick curves to show the continuous contributions to the 
AGB signal from sources at z greater than (0.5, 1.0, 1.6). 
We see that even if the popcorn and shot components can 
be identified, the continuous part of the AGB could mask 
a primordial GW background signal in the most sensitive 
frequency regime, around 20 - 30 Hz, of ET-B. As rates in- 
crease, so to does the continuous proportion of the AGB 
signal. 



9 CONCLUSIONS 

In this paper we have assessed the GW detection prospects 
for the background signals associated with SL-GRBs, as- 
suming that the central engines of a significant proportion of 
these bursts are provided by newly born magnetars. We have 
considered two plausible GW single-source emission mech- 
anisms: a) the deformation-induced GW emission from a 
newly born magnetar b) the onset of a secular bar-mode in- 
stability. The latter mechanism would correspond well with 
the long lived shallow plateau observed in the X-ray after- 
glows of many GRBs. 

We have calculated the GW background spectra of 
each of the mechanisms by employing appropriate models 
for each. We account for GRBs preference towards low- 




10 10 
frequency (Hz) 



Figure 11. As for Figure [8] but with thick lines showing only 
the continuous contribution (DC ^ 10) from sources (2 ^ 0.5, 
2 1.0, 2 > 1.6) for the rate estimates (ru, rp, r^). 



metallicity environments by using a source rate history 
model that allows for cosmic metallicity evolution. We as- 
sume both a low metallicity cutoff defined by e = 0.1 and a 
more relaxed cutoff, e = 0.5. 

Our results for the deformation-induced GW emission 
from a newly bor n magnetars are more pess imistic than 
those presented bv lRegimbau fc Mandid (|2008h . This is due 
to the fact that whilst they considered the emission from the 
population of magnetars assuming they represented 10% of 
the population of newborn neutron stars (a rate of around 
(3 — 15) x 10 3 Gpc _3 yr _1 ), we consider only magnetars as- 
sociated with SL-GRBs. 

For an AGB from triaxial emissions in newly formed 
magnetars associated with SL-GRBs, for an upper limit 
case in which emission is purely from GW emission, rates 
of (52, 137) Gpc~ 3 yr~ x will be required for detection at a 
SNR of 3 within 1 year by (ET-B, ET-C). For a AGB 
resulting from from toroidal fields, an upper rate ru = 
1800 Gpc _3 yr _1 would produce a SNR of 1.8 after 3 years 
of integration by ET-B. We find however that rates above 
200Gpc" 3 yr _1 would enable ET-B to explore the p B - B 
parameter space of the magnetar population considered in 
this study. 

In terms of detectability, we find that an AGB resulting 
from the onset of a secular instability is a more optimistic 
scenario. We note however, that this is highly dependent on 
rate of occurrence of thi s instability. Using the sin gle-source 
GW emission model of ICorsi fc Meszarosl (|2009al ). we find 
that event rates of (48, 80) Gpc~ 3 yr~ 1 are sufficient to pro- 
duce a detectable signal for ET (ET-B, ET-C) with SNR 
of 3 for 1 yr of observation. For ALIGO, detection within 3 
years would require the upper limit rate estimate ru, com- 
bined with a pair of optimally orientated and co-located 
detectors. We note that observations of a larger number 
of SL-GRBs (e.g. by future satellites like Janus or EXIST; 
IStamatikos et al.ll2009l ; llmerito et aLll2008l ). will help in re- 
ducing the uncertainties on their local rate estimates, thus 
clarifying the prospects of detectability of an associated GW 
background. 
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We find that this signal could potentially mask a pri- 
mordial GW background signal. Analysis of the DC for 
the background signal from secularly unstable magnetars 
showed that even at a conservative rate estimate tl, a sig- 
nificant proportion of this signal would be continuous. As 
highlighted in Fig. 1111 this would occur in the most sen- 
sitivity bandwidth of ET-B. Depending on the rate of oc- 
currence, both mechanisms could produce GW backgrounds 
that could mask a primordial GW background signal of or- 
der fis(/) ~ 5 x 10~ n in the frequency regime 10-50 Hz 
of ET-B and ET-C. This would pose a particular problem 
for the former, as it is the most sensitive bandwidth of this 
detector. 

The AGB may form a composite signal with an 
AGB from NS/NS inspirals, which is expected across 
a bandwidth of 10 - 800 Hz w ith increasing Q. b 
l|Regimbau fc de Freitas P achcco 2006). As the latter back- 
ground will peak at ~ 1000 Hz, detecting the higher fre- 
quency component may enable the two AGBs to be disen- 
tangled. 

To calculate SNRs we have applied equation (|22[) . which 
assumes the background signals can be detected through 
cross-correlation. We have chosen to adopt this convention 
for easy comparison with other astrophysical GW back- 
ground estimations, many of which will contain popcorn 
or shot noise components. In practice, to detect the non- 
continuous components of a GW background some other 
strategy will be required. Given the long duration and quasi- 
periodicity of the signal from secular instabilities, ET could 
detect a significant number of the shot noise events through 
matched filtering. This provides an additional means to in- 
terrogate the higher energy, shot noise and popcorn compo- 
nents of the background signal which result from the rarer 
nearby events. A statistical procedure, such as the "probabil- 
ity event horizon" technique, which extracts the observation 
time dependence fro m a population of cosmo logical tran- 
sients could be used l|Coward fc Burmanll2005T ). This tech- 
nique, which has been used to place constraints on the rate 
density of source populations, could use single detections 
from the shot noise component to interrogat e the tempo- 
ral dim e nsion of the GW background sig nal l|Howell et al.l 
2007tfi [C owar dl l2008l ; iHowell et al.ll2P/ITj i. 

Although the AGB signals discussed here remain spec- 
ulative, and the estimated rates are still highly uncertain, 
detection of the GW emission mechanisms associated with 
SL-GRB events could yield important payoffs. This could 
be possible should locally observed SL-GRBs produce the 
necessary triggers for multi-messenger observations. Cou- 
pled with a single-source detection, an AGB signal could 
provide constraints on the high-z evolution of the highly 
flux-limited SL-GRB population and would be a valuable 
probe of source rate and metallicity evolution. 



APPENDIX A: THE EFFECT OF COSMIC 
METALLICITY ON A GW BACKGROUND 
SIGNAL 

To illustrate the effect of metallicity dependence on calcula- 
tions of a GW background signal, in Figure lATl we reproduce 
the curves of Figs [6] and [8] for the rate rp. We add thick lines 
showing a metallicity independent (e = 1.0) source rate evo- 
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Figure Al. To illustrate the effect of metallicity dependence on 
the GW background signal we reproduce the curves of Figs 6 and 
[8] for rate rp. We add thick lines showing a metallicity indepen- 
dent source rate evolution assuming e = 1.0. We see that including 
a metallicity dependence shifts the background spectrum slightly 
to lower frequency. 

lution. We see that including a metallicity cutoff of e — 0.1, 
gives a higher signal at lower frequencies. This is a result of 
the greater contribution from high- z sources illustrated in 
Fig. [2] As indicated by the results in tables 3-6, this effect 
increases the SNR estimates. 
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